hep-th/0305022 



Classically integrable field theories with defects 



P. Bowcock a , E. Corrigan b and C. Zambon 

a Department of Mathematical Sciences 
University of Durham 
Durham DEI 3LE, U.K. 



b c 



Department of Mathematics 
University of York 
York YO10 5DD, U.K. 



ABSTRACT 



Some ideas and remarks are presented concerning a possible Lagrangian approach to 
the study of internal boundary conditions relating integrable fields at the junction 
of two domains. The main example given in the article concerns single real scalar 
fields in each domain and it is found that these may be free, of Liouville type, or of 
sinh-Gordon type. 
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1 Introduction 



Recently, there has been some interest in the study of integrable classical or quantum field 
theories restricted to a half-line, or interval, by imposing integrable boundary conditions, see for 
example [TJ El El IH El El Ej • The simplest situation, which is also the best understood, contains 
a real self-interacting scalar field with either a periodic (cos), or non-periodic (cosh) potential. 
The sinh-Gordon model can be restricted to the left half-line — oo < x < 0, without losing 
integrability, by imposing the boundary condition 



where m and (5 are the bulk mass scale and coupling constant, respectively, and e and e x are 
two additional parameters [21 E]- This set of boundary conditions generally breaks the reflection 
symmetry <ft — ► —<p of the model although the symmetry is explicitly preserved when e — £ i = £ - 
The restriction of the sinh-Gordon model to a half-line is a considerable complication, and 
renders the model more interesting than it appears to be in the bulk. This is because there will 
in general be additional states in the spectrum associated with the boundary, together with a 
set of reflection factors compatible with the bulk S-matrix (see El El EI)- The weak-strong 
coupling duality enjoyed by the bulk theory emerges in a new light [§1 ITU1 ITT]. 

In this article a slightly different situation is explored. There is no reason in principle why the 
point x = should not be an internal boundary linking a field theory in the region x < with 
a (possibly different) field theory in the region x > 0. The quantum version of this set up has 
been examined before and imposing the requirements of integrability was found to be highly 
restrictive. This sort of investigation was pioneered by Delfino, Mussardo and Simonetti some 
years ago J2], and there has also been some recent interest [T3*l 114j. However, the objective of 
this article is to explore a Lagrangian version of this question and derive the conditions linking 
the two field theories at their common boundary. This situation does not appear to have been 
discussed previously, although the results turn out to be interesting and reminiscent of some 
earlier work by Tarasov [To] . 

Internal boundary conditions will be referred to as 'defect' conditions. 

Integrability in the bulk sinh-Gordon model requires the existence of conserved quantities la- 
belled by odd spins s = ±1, ±3, . . . , and some of these should survive even in the presence of 
boundary conditions. Since boundary conditions typically violate translation invariance, it is 
expected that the 'momentum-like' combinations of conserved quantities will not be preserved. 
However, the 'energy-like' combinations, or some subset of them, might remain conserved, at 
least when suitably modified (see [2] for the paradigm). As was the case for the theory re- 
stricted to a half-line, the spin three charge already supplies the most general restrictions on the 
boundary condition. The Lax pair approach developed in can be adapted to this new context 
and used to re-derive the boundary conditions, thereby demonstrating that the preservation of 
higher spin energy-like charges imposes no further restrictions on the boundary conditions. This 
will be discussed below in section 4. 
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The starting point for the discussion is the Lagrangian density for the pair of real scalar fields 

01, 02: 



C = 0(-^)fI(a0O 2 -m))+^)U(90 2 ) 2 -W 2 )) (1-2) 

+ S { X ) Q Ol<%02 - 02<9 t 0l) - i3(0i, 2 ; 

in which the bulk potentials V\, V 2 depend only the fields 0i, 02, respectively, while the boundary 
potential B depends on the values of both fields at the boundary x = 0. The part of the boundary 
term depending on the time derivatives of the fields is not the most general possibility. However, 
although excluding terms of higher order in time derivatives, it is sufficiently general for present 
purposes. The field equations and associated boundary conditions are 

(1.3) 
(1.4) 

d x <pi - <9 4 2 = -^r, x = (1.5) 

>.r<>-2 ~ (/ /"J = • •' (1-6) 
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2 Consequences of the spin three conservation law 

For a single real scalar field in the bulk, the spin three densities satisfy 

d T r ±4 = d ± e T2 , (2.1) 



where 



T ±4 = A 2 (<9 ± 0) 4 + («9 2 ^ 



>± 



e ±2 = - 2 (d^Y w , (2.2) 

and (|2.1j) requires 

<9 3 ^ ,, 2 <9K . , 

W = AX W (2 ' 3) 

Thus, the only possibilities are a free massive field (A = 0, V — m 2 2 /2), a free massless field 
(V = 0, A ^ 0), or 
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with A, B being arbitrary constants. 

With two fields participating in different regions, the energy-like conserved quantities will be 
given by the following expressions: 



(2) 

-s+1 



+B S , (2.4) 
for a suitable boundary functional of 0, B s . The latter will be determined by requiring 
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(2.5) 



For the energy itself B\ = £>. For other values of s the argument is the familiar one from [2j, 
in the sense that the existence of B s and making use of (jOj) places severe constraints on the 
boundary potential B. Thus, for s = 3, and after some algebra: 
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(2.6) 
(2.7) 

(2.8) 
(2.9) 



There are several possible solutions to these constraints. The typical one, assuming neither of 
the fields is free and massive in its bulk domain, requires Ai = A 2 = A ^ 0. Consequently 
(ignoring an overall additive constant), 



B = ae x ^ 1+M + be x{<t>1 '^ ] + ce - ^ 1- ^ + de~ AW>1+<fc) , 

where a, b, c, d are constants, and the bulk potentials are given by 

V 1 = A x e 2X ^ + B ie ~ 2X ^ 
V 2 = A 2 e 2X<p2 + B 2 e~ 2X ^ 2 , 



(2.10) 



(2.11) 



with Ai = 2X 2 ab, A 2 = 2X 2 ac, B\ = 2X 2 cd, B 2 = 2X 2 bd. Notice that this case allows one of the 
bulk fields to be massless and free but the other need not necessarily be (for example, taking 
c = d = leads to a free massless field in x > 0, with a Liouville field in x < 0). 
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The alternative is that both fields are free and massive, so that Ai = A2 = 0. In that case, the 
conditions on the boundary and bulk potentials require the two masses to be the same, with a 
boundary potential of the general quadratic form: 

B = afil + &0102 + 

where a, b, c are constants. 

If one of the bulk fields is free and massless (say <pi), the other field may either also be free and 
massless, in which case the boundary potential has the form 

B = ae (01±02) +&- (01±02) , 

or the other field may be Liouville, in which case the boundary term has the form 

B = ™WV2 ( ae ^/V2 + I e -/WVA (2 12 ) 

P \ a J 

Finally, both fields could be of Liouville type. For example, choosing d — in (j2.1()j) leads to 
B\ — B 2 — 0, A\ 7^ 0, A 2 7^ and both potentials are Liouville potentials. 

Returning to the general case, the fields may be shifted by a constant in each bulk domain 
so that in ()2.11|) A\ = B\ and A 2 = B 2 , the latter in turn implying c = ±6, d = ±a. It is 
convenient to choose A = (3/ a/2, in order to agree with standard conventions for the sinh-Gordon 
model, and then to let a = ma/jS 2 , b = m/j3 2 cr. With those choices, the bulk and boundary 
potentials are: 

2 

771 ( V2(3<t> 2 , „-V2/3fc 



r 2 : — ( < 
B=^- (eP&i+fo)/^ ± e -P('t>i+<h)/V2' s j + ^(4>i-4>2)/V2 ± e -P(4>i-4>2)/>/2j ^.13) 

where, in all of these the '±' signs are strictly correlated. (In the sinh-Gordon model the relative 
signs cannot be adjusted by a real shift of one of the fields). There is a single free parameter a 
in the defect condition, which is perhaps puzzling since the half-line boundary condition ()2.1()j) 
allows two free parameters. 

Notice, the model might be restricted to a half-line if <p 2 say, were to be set to a constant value. 
Then the boundary condition satisfied by <pi would be of the general type (jl.lj) . However, 
typically, the constant value of <p 2 would not satisfy the equation of motion in 1 > 0. On the 
other hand, even if the equation of motion were to be satisfied with a constant <fi 2 in x > (ie 
for real (3, <fi 2 = 0), the boundary condition for 2 would generally not be satisfied at x = 0. 

Notice too, in none of these cases is there any reason why (f>i = (f> 2 at x = 0. Given the usual 
definition of 'topological' charge: 

r0 poo 

Q= dxd x (px + / dxd x (j) 2 = ^ 2 |oo - <pi\-oo + <f>i\o - 021c (2.14) 
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it is clear the difference 0i — 02 measures the strength of the defect at x = 0. In the sine-Gordon 
model, where similar considerations would apply, this would appear to indicate that topological 
charge need not be preserved, and hence that a defect need not necessarily preserve soliton 
number. 

It is worth noting that if eqs ljl.Hj) were satisfied simultaneously in the bulk then the 'defect' 
conditions with the choice of B given by ()2.10|) would become a Backlund transformation 16 
relating the two fields 0i and 02- Indeed, a sufficient condition for the Backlund transformation 
to work in the bulk would be: 

d 2 B d 2 B ( dB\ 2 ( dB 



% 90* ' 



204 - V 2 



Clearly both of these are satisfied in all the cases mentioned above. In the present setup, the 
'Backlund transformation' at x = represents the boundary between two domains. This sheds 
an interesting new light on the Backlund transformation itself. 

Generalising the idea, any number of defects may be represented similarly at domain boundaries 
x — Xi, x 2 , ■ ■ ■ , and at each boundary the defect conditions ought to retain the same form, albeit 
with different free parameters &i at each. In the bulk, the Backlund transformation between 
two solutions of the sine-Gordon equation generally changes soliton number (typically adding 
or subtracting a soliton), which appears to corroborate the suggestion above that a defect could 
allow a change of topological charge. In addition, different domains may contain fields of different 
character provided they are compatible with the boundary condition. 

One interesting further point. The canonical momentum density (which is not expected to be 
preserved because of the loss of translation invariance) is given by: 

/•O poo 
P = dx ^01^01 + / dx 8 t (j) 2 d x (j) 2 . 







Although P is not conserved, using the defect conditions (jl.Hj) it is not difficult to derive the 
following: 

dP ( . ± dB _ 1 dB 1 / dB . , , , , . , . . 

-dt<t>2^ r -d t <l) 1 — + -[ — \ - — -V 1 + V 2 \ . (2.15) 




dt V 90i "902 2 
The right hand side of ()2.15|) is a total time derivative provided that at x 
dB\ 2 ( dB\ 2 , d 2 B d 2 B 

These conditions are precisely satisfied by the boundary term indicated in ()2.10|) (and indeed co- 
incide with the conditions in the bulk mentioned earlier for a working Backlund transformation). 
In other words, there exists a functional of 0i, 2 , call it Vb, so that P ' + Vb is conserved. There 
appears to be a 'total' momentum which is preserved containing bulk and defect contributions. 
Thus, the fields can exchange both energy and momentum with the defect despite the lack of 
translation invariance. Clearly, there is a generalisation of this idea to a collection of defects 
situated at 
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3 The absence of reflection by defects 



Consider the consequences of a linearised version of (j2.13j) by setting 



e 



where R and T are reflection and transmission coefficients, respectively. (Strictly speaking the 
fields are the real parts of these expressions; however, in the linearised situation this is immate- 
rial.) Imposing the defect conditions, it is convenient to set k = 2m sinh uj = 2m cosh 9, and 
a = e p , to discover 

R{k)= 0, rw = _2ic°s^-(.-l/.) = _.sinh(^-f) 
KJ ' yJ 2ismh9-(a + l/a) sinh (2±£ + f ) V ; 

This is surprising, given the remarks about Backlund transformations, and should be compared 
with the results reported in [12] ■ There, the emphasis was different. The equations expressing 
the compatibility of reflection and transmission with the bulk factorisable S-matrix of a general 
model was found to be highly constraining and required the bulk S-matrix to satisfy S 2 = 1. 
In the context introduced here the question would be to find all the defect transmission factors 
compatible with the sinh-Gordon S-matrix, given that most probably there can be no reflection. 

To investigate what happens to a sine-Gordon soliton it is convenient to set (3 = l/y/2, m — 1/2 
and to write the bulk equations and defect condition as follows: 

x<0: <9 2 0i = -sin0!, (3.2) 
x > : <9 2 02 = -sin0 2 , (3.3) 

x = (J : o x <Pi — o t (f)2 = —0" sin | | sin 



2 J a 

o 1 rs, . f(pl+(p2\ 1 . f(j>l-(j>2\ f0 

o x (p2 - o t (pi = asm I J sin I 1. (3.4) 

Then, a single soliton solution in the two regions has the form (see for example [T7] ) 

e^ /2 = ^fk E a = C a e^\ al-f3l = l, a = 1,2, (3.5) 
1 + iE a 

where C a is real. In order to be able to satisfy the conditions ()3.4)1 the time dependence must 
match in the two domains (/?i = $2) and the constants C\, C2 are related by 

^ - (£f ) O., (3-6) 

where as before it is convenient to let a± = 02 = cosh 9 and /?x = P2 — sinh 9. Thus, the effect 
of the defect is to delay or advance the soliton as it passes through. One curious feature is that 
the defect can absorb or emit a soliton but only at a special value of rapidity. This is most 
easily seen by examining (|3.fi|) and noting that C2 vanishes for a < and e e = \a\, and C2 is 
infinite for a > and e e = a. In either case, the implication is that 02 = and a soliton with 
this special rapidity, approaching the defect from the region x < 0, will be absorbed by it. 
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4 Defect Lax pairs 



In this section the intention is to give an outline of the kind of approach one might adopt to set 
up Lax pairs in the presence of defects. 

To construct Lax pairs along the lines suggested in |7] it is necessary to separate slightly the 
boundary conditions in the two regions x < and x > 0, imposing the <pi boundary condition 
at x — a, and the 4>2 boundary condition at x = b > a, and to assume both fields are defined 
in the 'overlap' region a < x < b. Since the same framework applies to all the Toda, or affine 
Toda field theories this section will be quite general. Thus, with the same choices of coupling 
and mass scale as in the last section, the defect Lax pairs for models based on simply-laced root 
data are: 
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where for p = 1,2, 
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(4.2) 



Here H are the generators in the Cartan subalgebra of the semi-simple Lie algebra whose simple 
roots are «j, i = 1, . . . ,r, and E± a . are the generators corresponding to the simple roots or their 
negatives. If the theory is affine then the lowest root a = — n^ati is appended to the set of 
simple roots. In either case, affine or non-affine, the two expressions (14.2)1 are easily checked to 
be a Lax pair (for more details about this, and further references, see [TH]). 

To ensure the Lax pair defined by (J4.1)) really corresponds to a zero curvature in the overlap of 
the two regions there should exist a group element /C having the property 



d t JC = JCa ( o\t,b) -a { o\t,a)IC. 



(4.3) 



Setting 
with d t K 



fo-H/2 £ e <£vH/2 



has the effect of removing the time derivatives from the defect term in ([4.1)1 . 
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leading to 

0<pi 0<p 2 

(-\ e a ^ + ^> 2 [iC, E ai ] + i (e a ^- M/2 lCE„ ai - e M^i-^)/2 E^iC)) (4.4) 

i ' 

which is effectively an equation for both B and IC, see j7J. However, the structure of ()4.4|) is not 
the quite the same as that encountered previously. Nevertheless, a perturbative solution can be 
sought of the form 

K = 1 + _i + _f + ..., 

the 0(A) terms are identically satisfied, and the other terms lead to the following set of expres- 
sions 



0(1): (f + © H = -E^*^[*..« 

dB dB 
0(1/X): ^■Hk 1 + k 1 H--^- = 



>1 Uip2 



(4.5) 

The first of these can be satisfied for an arbitrary Toda model provided k\ = PiE- ai and 

however, the 0(1/ X) equation does not appear to be compatible with all choices of simple roots. 
Indeed, most Toda models appear to be ruled out. For the simplest, based on the A\ root 
system, a complete expression for JC is: 

^ = ( l) + J (l i ' Po = P 1 = P> ( 4 - 6 ) 



and 

B{<j> x -<h) = ~ (e a ^-^ 2 + e^-^ 2 ) . (4.7) 



P 

Here, the conventions used are: 



ai = a = V2 = -a , E a = = E -a-> H = (1/V2) ■ 

In other words, it appears this style of Lax pair can really only work in the presence of a defect for 
the sinh-Gordon or Liouville models. This result is identical with results obtained by examining 
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the spin two conserved charges for the models based on the A n root systems for n > 2, although 
these will not be reported in detail here. 

It was pointed out many years ago that the possible integrable boundary conditions are very 
constrained for all Toda models apart from the sinh-Gordon model in the sense that, in most 
cases, only a discrete set of parameters may be introduced at a boundary [HE). It now appears 
that defects are still more strongly constrained and generally cannot exist in the Lagrangian 
form postulated in ()1.2|) . 
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